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We analyse the rare decay mode B s — > fi~ in the minimal supergravity scenario (mSUGRA). We 
find a strong correlation with the muon anomalous magnetic moment (g — 2) M . An interpretation 
of the recently measured excess in (g — 2)^ in terms of mSUGRA corrections implies a substantial 
supersymmetric enhancement of the branching ratio B(B S — * fi + fM~): if (g— 2) M exceeds the Standard 
Model prediction by 4 ■ 10~ 9 , B(B S — > fi + fi~) is larger by a factor of 10-100 than in the Standard 
Model and within reach of Run-II of the Tevatron. Thus an experimental search for B a — > fi + fi~ 
is a stringent test of the mSUGRA GUT scale boundary conditions. If the decay B a — > /i + fi~ is 
observed at Run-II of the Tevatron, then we predict the mass of the lightest supersymmetric Higgs 
boson to be less than 120 GeV. The decay B s — » /i + fi~ can also significantly probe the favoured 
parameter range in SO(10) SUSY GUT models. 



Supersymmetry (SUSY) is an attractive and widely 
studied extension of the Standard Model (SM). The min- 
imal supergravity model (mSUGRA) Q relates all su- 
persymmetric parameters to just 5 real quantities: the 
universal scalar and gaugino masses Mq and M 1 / 2 , the 
trilincar term Aq , the ratio tan /3 of the two Higgs vac- 
uum expectation values, and sgn//, where [i is the Higgs- 
ino mass parameter. The first three quantities are de- 
fined at a high, grand unified energy scale and the others 
at the electroweak scale. They are the boundary con- 
ditions for the renormalization group equations, which 
determine the physical parameters at our low scale. Pre- 
cision observables, which are affected by SUSY correc- 
tions through loop effects, play an important role in 
constraining the supersymmetric parameter space. The 
small number of parameters makes mSUGRA highly pre- 
dictive so it can be significantly tested by low energy 
precision measurements. In this letter we show that the 
decay B s — > /J + /i~ is a stringent test of the mSUGRA 
scenario, in particular when correlated with [g — 2) M . 

Recently the Brookhavcn National Laboratory (BNL) 
reported an excess of the muon anomalous magnetic mo- 
ment = (g - 2)^/2 over its SM value |g. The differ- 
ence 6a M =a e *V-al M = (43 ± 16) • 1(T 10 corresponds to 
a 2.6a deviation from the SM. An mSUGRA interpreta- 
tion of this anomaly implies /i > (in the sign convention 
with Mi/2 > an d equal signs of the diagonal elements 
of the chargino mass matrix) J3|. It further invites a 
large tan/3 > 10 Q|. The discrepancy in the case of a M 
is by itself not significant enough to justify the claim 
of new physics, especially since the calculation of a^ M 
involves two hadronic quantities: the hadronic contribu- 
tions to the photon self-energy, which must be obtained 
from other experiments, and the (smaller) light-by-light 
scattering contribution, which can only be estimated with 
hadronic models. A more conservative estimate of the 
latter would reduce the BNL anomaly to a 2a effect [||. 
Hence in order to resolve the possible ambiguity between 
mSUGRA and alternative explanations of Sa^ one ideally 
wishes to study other observables whose sensitivity to su- 



persymmetric loop corrections is correlated with Sa^. It 
is our purpose here to show the strong correlation be- 
tween B(B S — > /j + /i~) and Sa^ in mSUGRA. 

SUSY modifies B meson observables if tan/3 is large, 
because the b Yukawa coupling becomes sizable. Espe- 
cially sensitive are quantities with a b quark chirality flip 
like the branching ratios B[B -> A s7 ) and B{B -> £+£~). 
In mSUGRA the low energy value for the trilinear term 
A t is dominated by M 1/2 with A t < for M 1/2 > [pj . 
Then > implies that the charged-Higgs-top loop and 
the chargino-stop loop tend to cancel in B(B — > X s -f), so 
that the sensitivity to mSUGRA corrections is weakened. 
A further disadvantage of this decay mode is that it re- 
quires an experimental cut on the photon energy, which 
introduces some hadronic uncertainty. 

In [0 the possible impact of flavour-blind SUSY on 
other B physics observables, in particular those which en- 
ter the fit of the unitarity triangle, were studied and only 
small effects were found. This did not include the decay 
B s — ► In contrast to the observables in @], the 

branching ratio B(B S — > /j + /i~) grows like tan 6 /3 [|| |J, 
with a possible several orders of magnitude enhancement. 
We here go beyond this work to study B s — > in 
the mSUGRA model. Since B(B -> £+£~) oc mj, the 
branching ratio is largest for £ = t. Yet r's are hard to 
detect at hadron colliders, so that the prime experimen- 
tal focus is on the search for B s — > B factories 
running on the T(45) resonance produce no B s mesons. 
Leptonic branching ratios of Bd mesons are smaller by a 
factor of \Vtd/Vts\ 2 ^ 0.06. Since in B factories the boost 
of the Bd meson is known and the considered leptonic 
decay rates can be substantially enhanced over their SM 
values in SUSY, we encourage our colleagues at BaBar 
and BELLE to look for Bd — ► t + t~ decays, as well. 
From now on we restrict ourselves to the decay mode 
B s ^JJ-^J- 1 - 

In Ji ll the SUSY corrections to B{B S -> m + aO were 
calculated at the one-loop level. For large tan/3, higher 
order corrections can be large, eventually of order 1. 
In P] tan /3-enhanced supersymmetric QCD corrections 



1 



have been summed to all orders in perturbation the- 
ory. We have incorporated these dominant higher or- 
der corrections by replacing the b Yukawa coupling hb oc 
nib tan/3 with ft| ff = hb/(l+Amb), where Am& oc fi tan /3 
depends on the gluino and sbottom masses and can be 
found in ||. Amb is positive for \i > 0. The dominant 
contribution to B(B S — » /t + /U~) is proportional to hf 14 , 
so that the inclusion of Amb tempers the large-tan /3 be- 
haviour. 

The considered branching ratio can be expressed as 
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approximate formula is not accurate. A similar esti- 
mate of the supersymmctric contribution to yields 
(Sa^susY oc tan /3/ (M ) /A/^. (<Sa„) SUSY depends on 
slepton masses, which are less sensitive to M 1 / 2 than 
squark masses; they are dominated by M . We have en- 
coded the M dependence in the slowly varying function 
/(M ). Hence both B(B S -> |U + /x~) and (Sa^suSY 
grow with tan/3 and decrease with increasing M 1 / 2 . For 
this it is essential that we have made the assumption 
of the mSUGRA GUT scale boundary conditions. Thus 
within mSUGRA we expect a strong correlation between 
B{B S — > and (Sa^susV- if Sa^ ^ requires a 

SUSY explanation with large tan (3 then we would expect 
B(B S — + /t + /i~) to be strongly enhanced. If, however, the 
supersymmetric explanation of <5a M requires small M\i 2 
and a moderate value of tan (3 then we would expect only 
a moderate enhancement of B(B S — > /i + /i~). 



Here |Vt s | = 0.040 ± 0.002 is the relevant CKM ma- 
trix element and f Bs = (230 ± 30) McV @ is the 
B s decay constant. In (|l|) we have kept the depen- 
dence on the lepton mass m„, so that the generalisa- 
tion to B c i — > t + t~ is straightforward. The Wilson 
coefficients Cg, Cp and Ca, which contain the short- 
distance physics, are normalised as in O]. The co- 
efficients cs, cp and cio defined in 13] are related to 
ours by C,s = — 2cs sin 2 9yy, Ca = — 2cio sin 2 6w and 
C P = 2cpsm 2 e W - Within the SM, C s and C P are 
negligibly small and the NLO result for Ca can be 
well approximated by C A = 2.01(m t /167 GcV) 1 ' 55 Q. 
Here m* = rfitimt) is the top quark mass in the MS 
scheme. m t = 167 GeV corresponds to a pole mass 
of m t = 175 GeV. The SM prediction is given by 
B(B S -> /i+aO = (3.7 ± 1.2) • 10~ 9 , with the uncer- 
tainty (±25%) dominated by Jb b ■ This is also the main 
hadronic uncertainty in the SUSY calculation. 

During Run-I of the Tevatron, CDF determined E3] 



B{B S -> < 2.6 x 10" 



at 95% C.L. 
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The single event sensitivity of CDF at Run-IIa is es- 
timated to be 1.0 • 10~ 8 , for an integrated luminosity 
of 2 fb _1 @. Thus if mSUGRA corrections enhance 
B(B S — ► /i + /i~) to e.g. 5 • 10~ 7 , one will see 50 events 
in Run-IIa. Run-IIb may collect 10-20 fb _1 of integrated 
luminosity, which implies 250-500 events in this example. 

In SUSY, the dominant coefficients are Cs.p since 
they are proportional to tan 3 /3. We would like to un- 
derstand the effect of the restricted mSUGRA param- 
eters on Cs,p and thus on B(B S — > /i+/i~). In 
mSUGRA the low-energy values of both fj, and the 
squark masses are dominated by the (GUT scale) value 
of Mi/2 through the renormalization group equations. 
For not-too-large M ,M 1/2 < 500 GeV and A ~ 
GeV we can derive the approximate formula B(B S — > 
(i+fjL-) w 10- 6 tan 6 /3M 2 /2 GeV 4 /(M 2 /2 + Af 2 ) 3 . In the 
vicinity of the maximum (near M\/ 2 = 0.4 Mq) the 
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FIG. 1. (Sa^susY, versus B(B S — > n + fi~) for tan/3 (top) 
and Af 1/2 =450, M = 350, 4 = 0,fJ, > 0,m t = 175 
GeV. Shown also, the SM prediction, the present bound by 
CDF [|l3l, on B(B S — ► as well as the present lcr and 

We used /a. = 230 MeV. 

We now study these effects quantitatively. For this 
we use the full computation of Eq.(|l|) including the rc- 
summcd SUSY QCD corrections, and restricting our- 
selves to the mSUGRA parameters. In Fig.|l], we 
show the direct correlation between B(B S — > n + n~) 
and {5a^)susY for the fixed parameters: Mi/ 2 = 
450 GeV, Mo = 350 GcV,A = 0,/x > and m t = 175 
GeV. On the upper edge we show the tan/3 dependence. 
We restrict ourselves to tan /3 < 58 in order to guarantee 
radiative electroweak symmetry breaking (REWSB). We 
have included the SM prediction and the CDF bound 
from Eq.(^). The solid (dashed) curve represents the 
B(B S — > /i + /j~) result with (without) rcsummation of 
the tan /3-enhanccd SUSY-QCD corrections. In this ex- 
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ample, the resummation suppresses B(B S — > fi + fi~) by 
75% for tan/3 > 50. In order for mSUGRA to account 
for <5a M within la of the current BNL measurement at 
this parameter point, we see that we need a large value 
of tan/3>50. Due to the strong correlation within 
mSUGRA we then predict B(B S -> m + M~) £ 5 • 10~ 8 , 
which is observable by CDF at Run II. 
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FIG. 2. Contours of 6(£? s — > /i + /i _ ) (solid) and [Sa^susY 
(in units 10~ 10 ) (dashed) in the Mi/2-tan/3 plane. The light- 
est neutral CP-even Higgs mass is shown as well (dot-dashed). 
The shaded regions are excluded, as described in the text. 
The mSUGRA parameters are given at the top. 
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FIG. 3. Contour plots of the B(B S — » n + n~) (solid) and 
on (5a^)susY (dashed) in the (Mo, Mi/ 2 )-plane for mSUGRA 
parameter values as shown. The shaded regions are excluded 
as described in the text. Contours of the light Higgs boson 
mass (dot-dash line) are also shown. 

As we discussed above, we expect B(B S — > /i + /i~) to 
dominantly depend on the mSUGRA parameters Mj/2 
and tan/3. In Fig.|| we show the B(B S — ► (solid) 



and the (Sa^susY (dashed) contours in this plane. We 
have fixed: M = 300 GeV, Aq = 0, [i > and 
m t = 175 GeV. The 2er contours for Sa^ (11,75) are 
explicitly given. The left vertical shaded region is the- 
oretically excluded since it does not allow for REWSB 
or violates the LEP chargino bound. The upper right 
triangular shaded region is excluded, since the LSP is 
not neutral. If as expected, CDF can probe down to 
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") > 2-10" 7 at RUN Ha, this corresponds to 



a sensitivity of (Mi/ 2 , tan /3) ranging from (160 GeV, 47) 
to (450 GeV, 57). The qualitative discussion of before 
is now nicely reproduced. B(B S — > decreases 
with increasing M1/2 and rapidly increases with tan/3. 
Fig.^ also nicely shows the cross-correlation between 
{Sa^suSY and B(B S -> If both B(B S -> ^ + M~) 

and (5a fj,) susy are found in disagreement with the SM 
and are measured with a 50% and 20% accuracy, respec- 
tively, then for given Mo, this fixes tan/3 to better than 
20% and M 1/2 to better than 30%. 

It is conventional to discuss mSUGRA physics in 
the (M1/2, Mo)-plane. In Fig.|| we show the contours 
of B(B S — > (solid) and (Sa^susY (dashed) in 

this plane, for tan/3 = 50, Aq = 0, fj, > and 
m t = 175 GeV. Again we include the CDF bound 
Eq.(|[) and the Higgs mass contours. The left shaded 
region is excluded through the requirement of REWSB 
or the chargino bound. The lower right shaded re- 
gion is excluded through the requirement of a neutral 
LSP. A sensitivity of B{B S -> H+fi') > 2 ■ 10~ 7 at CDF 
now corresponds to a sensitivity of M 1 / 2 iS 280 GeV and 
M < 400 GeV, respectively. 

While CDF is not able to see squark masses directly up 
to 0.7 TeV (corresponding to M 1/2 = M ~ 300 GeV,), 
it will nevertheless be able to prepare the ground for 
LHC by observing the B s — » mode. Even bet- 

ter, after lOffcT 1 CDF will probe M 1/2 < 450 GeV and 
M < 600 GeV (for tan (3 = 50) which in mSUGRA cor- 
responds to masses for the heaviest superpartners of 1 
TeV. We conclude the discussion of Fig.|| with the pre- 
diction of the light Higgs boson mass Mh (dot-dashed 
line) for tan/3 = 50 in the mSUGRA scenario 15 1. Any 
measurement of B(B S — > by itself implies a use- 

ful upper bound on M^. The simultaneous information 
of B(B S —> /i + /i~) and <5a M fixes M/, in most regions of 
the (M 1 / 2 , M )-plane. A Higgs mass around 115.6 GeV 
results in 10~ 8 < B(B S -> n+fx') < 3 ■ 10~ 7 which would 
most likely be measured before the Higgs boson is dis- 
covered. 

In Figs. 1-3 we have chosen Ao = 0. A non-zero Aq 
changes the value of A t at low energies. This parame- 
ter plays a crucial role for the GIM cancellations among 
the contributions of different squarks to B(B S —> /i + /i~). 
Changing Ao to —500 GeV in the scenario of Fig. 1 en- 
hances B(B S — > /j + /i~) by up to a factor of 6 com- 
pared to the case with Aq = 0. For Aq = +500 GeV 
B(B S — > /i + /i~) is slightly decreased. 
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In our figures we have omitted further constraints on 
the mSUGRA parameter space, in order to clearly show 
the correlation between B(B S — > and {Sa^susY ■ 

The most significant further constraint comes from the 
measurement of B(B — > X s j) whose prediction is 
less certain in the large tan/? region jl7],D ■ If we take 
the conservative approach of jlS ] , then we can exclude 
values of M 1/2 < 250 GeV in Fig. | for tan [3 > 25. In the 
scenario of Fig.|| this implies B(B S — > < 5 ' 10 -7 . 

For a discussion of the constraints from supersymmetric 
dark matter see for example and references therein. 

The large values of tan /3 we have been considering 
are theoretically well motivated within SUSY SO(10) 
Yukawa unification. There a narrow parameter region can 
explain the observed Sa^ while still being consistent with 
the constraint from b — > sy 1 19 2^ . This is not within the 
context of mSUGRA. However, in this parameter region 
both fi and Mx/ 2 arc light, while the CP-odd Higgs boson 
mass is less than 300 GeV, and tan/3 w 50. Therefore 
we expect B(B S — > to be strongly enhanced. As 

an example we determine B(B S — > /i + /i~) for the best fit 
points found in @: A/ A = 110 GeV, m Sl < 250 GeV, 
| A t | > 1 TeV, m~ t < 1 TeV and tan/3 ~ 50. Within the 
hadronic uncertainties B(B S — > > 10~ 5 which is 

already excluded by CDF @. Thus the SO(10) models 
should be reconsidered in the light of B(B S — ■> /i + /i~). 
Turning it around, if an SO(10) GUT model is the cor- 
rect description of nature then the decay B s — > 
must be just around the corner. 

In conclusion, we have found a striking correlation be- 
tween the muon anomalous magnetic moment a M and the 
branching ratio B(B S — > fi + fi~) in mSUGRA scenarios. 
If the reported excess in M is caused by mSUGRA 
corrections with large tan/3, one faces more than an or- 
der of magnitude enhancement of B(B S — » over its 
SM value. This is within reach of Run-II of the Tevatron. 
The combined measurements significantly constrain the 
mSUGRA parameters, allowing a determination of tan (3 
and Mi/2- A measurement oiB(B s — ► will further 
constrain the mass of the lightest Higgs bosons. SO(10) 
SUSY explanations of the measured a M are barely com- 
patible with the present upper bound on B(B S — > 
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